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Abstract 


Diapoma pampeana was recently described to occur in the upper Negro basin in Uruguay and Brazil. An isolated population ten- 
tatively identified as D. pampeana from the Pando stream, a perturbed coastal drainage in Uruguay, is studied and compared under 
the light of morphological and molecular data to test if there is evidence to consider it as a separate species. New geographical re- 
cords for the species are presented and included in the comparisons. The specimens analyzed were pooled into four groups: Pando, 
Santa Lucia, Middle Negro and Upper Negro. We analyzed 32 morphological characters using statistical procedures and recovered 
a COI-based phylogeny of different populations of D. pampeana to test if they may represent different species. Size-corrected PCA 
revealed that the Pando and Upper Negro groups are greatly diverging in both morphometric and meristic data along PC1 (mainly 
by the snout to dorsal-fin origin, dorsal to adipose-fin origins, number of longitudinal scales and predorsal scales). This deviating 
pattern was also obtained in a cluster analysis. The Santa Lucia and Middle Negro groups were found to be intermediate morpho- 
types. In contrast, molecular analyses revealed that the Pando and Upper Negro specimens resemble genetically and, thus, are 
placed together in the Neighbor-joining and Bayesian topologies, as part of a monophyletic Diapoma. We proposed that the Pando 
population, despite its deviating morphology observed, can be classified as D. pampeana. Therefore, this population constitutes a 
remarkable example of an isolated population that is morphologically divergent but genetically similar to the geographically most 
distant conspecific population. 
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Introduction small-sized species (no more than 100 mm SL) (Thomaz et 
al. 2015; Mirande 2019; Ferreira et al. 2021; Ito et al. 2022; 
The Neotropical fish genus Diapoma Cope, 1894 isamem- Fricke et al. 2023). Diapoma is recognized as a monophy- 
ber of the tribe Diapomini, which is one of the largest mono- _letic group based on molecular data and combined evidence 
phyletic groups within the Stevardiinae with ~ 135 relatively (Thomaz et al. 2015; Mirande 2019; Ito et al. 2022). 
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To date, this genus includes sixteen valid species that 
are distributed along different river drainages in Argen- 
tina, Brazil, Paraguay and Uruguay, mainly within the 
Rio de la Plata basin. Five species are known from the 
Parana-Paraguay system, D. guarani (Mahnert & Géry, 
1987) (several streams flowing into the Parana basin in 
the border region between Alto Parana, Paraguay and Mi- 
siones, Argentina), D. obi (Casciotta, Almiron, Pialek & 
Rican, 2012) (some tributaries from the Paranay-Guazu 
drainage and the Moreno stream in Misiones, Argentina), 
D. nandi Vanegas-Rios, Azpelicueta & Malabarba, 2018 
(the Piray-Mini stream in Misiones, Argentina), the re- 
cently described D. potamohadros Ito, Carvalho, Pava- 
nelli, Vanegas-Rios & Malabarba, 2022 (endemic from 
the Rio Iguazu basin in Argentina and Brazil) and D. ter- 
ofali (Géry, 1964) (the Rio Lujan and other streams flow- 
ing into the Rio de la Plata basin in Buenos Aires, Argen- 
tina) (Mahnert and Géry 1987; Menezes and Weitzman 
2011; Casciotta et al. 2012; Vanegas-Rios et al. 2018; Ito 
et al. 2022). The Rio Uruguay basin is the water body that 
possesses the greater number of Diapoma species regis- 
tered so far (seven spp.), D. alegretense (Malabarba & 
Weitzman, 2003) (the Rio Ibicui system in Rio Grande do 
Sul, Brazil, and Laguna Redonda in Artigas, Uruguay), 
D. guarani (Barragem Sanchuri in Rio Grande do Sul, 
Brazil), D. lepiclastum (Malabarba, Weitzman & Casci- 
otta, 2003) (from the Rio Pelotas and Rio Canoas to the 
Lageado Uniao stream in Brazil, and in the eastern region 
of Misiones, Argentina), D. pampeana Ito, Carvalho, 
Pavanelli, Vanegas-Rios & Malabarba, 2022 (the upper 
Rio Negro basin in Brazil and Uruguay), D. pyrrhopteryx 
Menezes & Weitzman, 2011 (the Rio Pelotas basin in Rio 
Grande do Sul, Brazil, and the Pepiri Iguazu basin in Mi- 
siones, Argentina), D. terofali (streams flowing into the 
Rio Uruguay system in Rio Grande do Sul, Brazil, and in 
Artigas and Cerro Largo, Uruguay), and D. uruguayense 
(Messner, 1962) (tributaries of the Rio Uruguay in the 
border region between Argentina, Uruguay, and Brazil, 
and from the headwaters of the Rio Negro) (Malabarba 
and Weitzman 2003; Zarucki et al. 2010; Menezes and 
Weitzman 2011; Thomaz et al. 2015; Almiron et al. 2016). 
Four species have been recorded for the Laguna dos Patos 
basin and the coastal drainages of south Brazil, D. dicrop- 
otamicum (Malabarba & Weitzman, 2003) (northern trib- 
utaries of the Rio Jacui from the Serra Geral formation), 
D. itaimbe (Malabarba & Weitzman, 2003) (Tramandai, 
Mampituba, Ararangua river basins, southern coast of 
Brazil), D. speculiferum Cope, 1894, D. thauma Menez- 
es & Weitzman, 2011 (tributaries of the Rio Jacui basin, 
Rio Grande do Sul, Brazil) and D. tipiaia (Malabarba & 
Weitzman, 2003). Finally, D. alburnum (Hensel, 1870) 
is the most widely distributed species occurring in the 
Uruguay (e.g. the Rio Queguay basin), Parana (the Rio 
Gualeguay basin), Laguna dos Patos (e.g. the Rio Jacui 
system) basins and coastal drainages of southern Brazil 
(Malabarba 1983; Malabarba and Weitzman 2003; Proto- 
gino and Miquelarena 2012; Paullier et al. 2019). 

Diapoma pampeana, recently described from the 
upper Rio Negro basin, reaches 35 mm SL and can be 
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differentiated from all its congeners by a combination 
of characters, mainly from body pigmentation (Ito et 
al. 2022), including the presence of: a narrow and con- 
spicuous black line along the horizontal septum, never 
forming a wide lateral stripe; a longitudinal black stripe 
extending posteriorly on the middle caudal-fin rays; and 
a small black blotch, restricted to the base of the middle 
caudal-fin rays. 

We found specimens that potentially could be identified 
as D. pampeana from the Pando stream, a coastal drain- 
age flowing into the Rio de la Plata estuary in Uruguay, 
based on the resemblance of the humeral mark, midlateral 
stripe, and caudal-fin pigmentation. The possible pres- 
ence of this species in the Pando stream caught our atten- 
tion because the preliminary morphometric data obtained 
were somewhat incongruent with the data reported in the 
description by Ito et al. (2022), and because no conge- 
ner has been recorded in this area so far (Malabarba and 
Weitzman 2003; Menezes and Weitzman 2011; Gurdek 
and Acufia-Plavan 2017). This small drainage, which 
is located in the middle of urban and agricultural areas, 
has been greatly modified and affected by anthropogenic 
factors such as pollution, industrial activities, and urban 
waste originated from anthropogenic actions (Echevar- 
ria et al. 2011; Achkar et al. 2012; Gutiérrez et al. 2015; 
Muniz et al. 2019) and is considered of great importance, 
among other reasons, because its sub-estuarine mouth 
plays a role in the breeding and nursery for grounds of fish 
(Defeo et al. 2009; Acufia et al. 2017; Muniz et al. 2019). 

There are several studied cases in which species pre- 
viously considered as distributed in the Rio Uruguay 
drainage and Atlantic river coastal drainages have been 
separated in two different species [e.g. Parapimelodus 
nigribarbis (Boulenger, 1889) vs. P. valenciennis (Lut- 
ken, 1874), see Lucena et al. (1992); Pimelodus pintado 
Azpelicueta, Lundberg & Loureiro, 2008 vs. P. maculatus 
Lacepeéde, 1803, see Azpelicueta et al. (2008); Bunoceph- 
alus erondinae Cardoso, 2010 vs. B. doriae Boulenger, 
1902, see Cardoso (2010); Pseudocorynopoma stanleyi 
Malabarba, Chuctaya, Hirschmann, Oliveira & Thomaz, 
2020 vs. P. doriae Perugia, 1891, see Malabarba et al. 
(2021)]. So, the present study aims to carry out a mor- 
phological (mainly morphometric and meristic data) 
and COlI-based comparison between the population of 
D. pampeana from the Rio Negro basin, a tributary of the 
lower Rio Uruguay, and the population tentatively identi- 
fied as D. pampeana from the Pando stream, that empties 
directly in the Atlantic Ocean. We expect to test if there is 
evidence to consider these isolated populations as sepa- 
rate species or to treat them as a single species, describing 
any intraspecific variation between them. Additionally, a 
recent examination of specimens of Diapoma at Museo 
Nacional de Historia Natural, Montevideo (MHNM), 
revealed two lots of individuals similar to D. pampea- 
na in the body shape and meristic data from the Yi and 
Santa Lucia river basins (as fixed in 10% formalin, DNA 
extraction was unavailable). Consequently, they were in- 
cluded in the morphological analyses to enhance the com- 
parisons and are also presented as new records. 


Zoosyst. Evol. 100 (1) 2024, 69-85 


Materials and methods 


Specimens were collected in the Pando stream (permis- 
sion No. 202/717/04, DINARA, Uruguay) between 2003 
and 2004. They were fixed in formalin 10% and preserved 
in alcohol 70%. Some of them, which were preserved 
originally in ethanol 96%, were rehydrated before being 
preserved and catalogued as the others. Additional stud- 
ied specimens of D. pampeana and comparative species 
are deposited in the following institutions: MACN-ict, 
MLP, MHNG, MHNM, UFRGS, and UNMDP (abbrevi- 
ations according to Sabaj 2020). 


Morphological analysis 


Measurements and other counts were taken following 
Fink and Weitzman (1974), with the modifications pre- 
sented by Ito et al. (2022). Twenty-two measurements 
were taken point to point with a digital caliper under a 
stereomicroscope and are expressed as percentages of 
standard length (SL) or head length (HL) for units of the 
head. Specimens were cleared and stained (c&s) follow- 
ing Taylor and Dyke (1985). The total number of ver- 
tebrae was counted in c&s specimens. Those counts in- 
cluded the first preural centrum plus the first ural centrum 
(PU1 + U1) counted as one element and all four vertebrae 
of the Weberian apparatus. 

The specimens from the Pando stream were compared 
with type specimens of D. pampeana from the Rio Negro 
basin under different statistical procedures. Additionally, 
the specimens presumably belonging to D. pampeana 
from the Yi and Santa Lucia river basins in Uruguay were 
only processed in the morphological analyses because 
they were not suitable for DNA extraction. To facilitate 
comparisons, a morphometric data matrix that included 
all of these specimens was pooled into groups (based on 
geographic drainages) as follows: Pando (n = 17), Santa 
Lucia (Canelon Grande, n=2), Middle Negro (Y1, n= 15), 
and Upper Negro (several localities, n = 35). Nearly all 
the specimens analyzed in all groups were adults, except 
for a few immature specimens that were also added in the 
comparisons (excluding the bone hooks, no morphomet- 
ric or meristic differences were observed between them 
and the respective adults). This dataset was analyzed 
using the “allometric vs. standard” procedure (Elliott et 
al. 1995), under which the allometric coefficients are cal- 
culated concerning a standard (reference, such as over- 
all length) measurement (each variable is regressed onto 
this after log-transformation). The size-corrected mor- 
phometric dataset was analyzed using a principal com- 
ponent analysis (PCA), based on the covariance matrix. 
For the PCA, the number of significant principal compo- 
nents (PCs) was decided by two criteria: the broken-stick 
model (Frontier 1976) and the scree plot method (Cattel 
1966). To compare the dissimilarity between the groups 
associated with the size-corrected morphometric data, a 
hierarchical cluster analysis was performed using Ward’s 
method (Ward 1963) and Euclidean distances, under 1000 


71 


bootstrap replicates. Missing values in measurements 
(e.g. some fin rays were broken in a few individuals) were 
imputed from predictor values obtained under maximum 
likelihood from the EM algorithm (Dempster et al. 1977; 
Pigott 2001) using 500 iterations. For the morphometric 
data, confidence intervals of 95% were calculated using 
9000 bootstraps. 

Meristic data showing slightly different patterns be- 
tween the Pando group and the other groups were analyzed 
using a PCA on the root-squared transformed values and 
the correlation matrix (Quinn and Keough 2002) (outli- 
ers that could not be reexamined in the specimens were 
omitted and the mean was imputed for missing data). To 
illustrate the distinctive patterns within the Pando group, 
we used Tukey box plots for the variables, which provide 
a clearer representation of the observed variability. When 
it comes to counts, both mean and mode values are re- 
ported, and they are separated by a slash. 

For those analyses, normality was tested using a Sha- 
piro—Wilk statistic (W) in each case (a < 0.05) and data 
were log-transformed when needed to better approximate 
to a multivariate normality. Statistical procedures were 
carried out in PAST 4.12 (Hammer et al. 2001), IBM 
SPSS Statistics 26.0 (IBM 2019), and GraphPad Prism 
9.4.1 (GraphPad Software, San Diego, CA, USA). 


Molecular analysis 


The mitochondrial cytochrome c oxidase subunit I gene 
(COI) was obtained from two specimens from the Pando 
stream. DNA extraction and polymerase chain reaction 
(PCR) were carried out following the standard COI proto- 
cols (Ivanova et al. 2006; Rosso et al. 2012), under differ- 
ent sets of primer cocktails for fishes (Ivanova et al. 2007). 

In total each amplification reaction produced a vol- 
ume of 12.375 uwL from 2 wl of DNA template, 6.25 wL 
of 10% trehalose, 2 uL of molecular biology grade wa- 
ter, 1.25 wL of 10x reaction buffer, 0.625 wL of MgCl2 
(50 uM), 0.0625 wL of dNTP (10 mM), 0.0625 wL of each 
primer (10 uM) and 0.0625 uL of Invitrogen’s Platinum 
Taq. polymerase (5 U wL-1). The amplification condi- 
tions consisted of 2 min at 95 °C, followed by 35 cycles 
at 94 °C for 30 s, at 52 °C for 40 s and at 72 °C for 1 min, 
and ended at 72 °C for 10 min. E-Gels (Invitrogen) were 
used to check the amplification success. The COI gene 
was sequenced in Macrogen (Korea) and IGEVET-UNLP 
(Argentina). Sequence chromatograms were edited using 
BioEdit 7.2.5 (Hall 1999). 

For comparative purposes, in addition to the newly 
generated sequences, 72 COI sequences were selected 
from representative specimens of the valid species of Dia- 
poma (except D. nandi) (Ito et al. 2022) and genera close- 
ly related to it (detailed later herein) that are available in 
GenBank and Barcode of Life database (BOLD, available 
at http://www.boldsystems.org) (accession numbers for 
all sequences analyzed are provided in Suppl. material 1). 
The COI sequences were aligned with MUSCLE (1000 
iterations) (Edgar 2004) and generated as a data matrix 
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partitioned by the first three codon positions in MEGA 
11.0.13 (Tamura et al. 2021). The COI dataset was upload- 
ed to Zenodo (https://do1.org/10.5281/zenodo.8361520). 

To analyze the phylogenetic placement of the specimens 
from the Pando stream through different methods, the COI 
data matrix was analyzed by the phylogenetic procedures 
and conditions described hereafter. Modeltest-NG (Darri- 
ba et al. 2019) was used for selecting the best-fit nucleo- 
tide substitution model available for each procedure (and 
computational package) based on the partitioned align- 
ment when necessary. The AIC and BIC statistical criteria 
were explored, but the latter was used to choose the best 
model among the candidate models. The neighbor-joining 
(NJ) tree (10000 bootstrap) was constructed based on the 
Tajima+Nei model with rates gamma-distributed as im- 
plemented in Mega. Bayesian analyses were conducted in 
MRBAYES 3.2.2 (Ronquist et al. 2012) using two runs, 
each with four Markov chains, which ran for 60 million 
generations (25% discarded as burn-in, sampling a tree 
every 3000 generations) under the models SYM+I (posi- 
tion 1), F81 (position 2) and HK Y+G (position 3). Tracer 
1.7.2 (Rambaut et al. 2018) was used to evaluate the re- 
sults of the MRBAYES analyses from each run (ESS val- 
ues and InL plots). The trees were visualized and prepared 
with FIGTREE 1.4.4 (Rambaut 2018). The CIPRES por- 
tal (Miller et al. 2010) was used to run the following com- 
putational programs: MODELTEST-NG and MRBAYES 
3.2.2. The outgroup was composed of species from Bry- 
conamericus Eigenmann, 1907, Hypobrycon Malabarba 
& Malabarba, 1994, Nantis Mirande, Aguilera & Azpeli- 
cueta, 2006, Odontostoechus Gomes, 1947, Piabarchus 
Myers, 1928, and Piabina Reinhardt, 1867, genera that 
have been considered in preceding studies (Ferreira et al. 
2011; Mirande 2019; Ito et al. 2022) to be closely related 
to Diapoma (trees were rooted in Piabina species when 
required). Additionally, to compare the interspecific and 
within-species variability, the uncorrected pairwise ge- 
netic distances (gamma distributed and pairwise deletion) 
were calculated in MEGA 11.0.13 (Tamura et al. 2021). 

To examine the potential variability associated with 
polymorphism between the Pando and Upper Negro 
specimens, despite the limited number of samples avail- 
able, the specimens of D. pampeana in the COI data 
matrix were pooled into the two respective groups using 
DNASP 6.12.03 software (Rozas et al. 2017) to calculate 
the polymorphism sites, nucleotide diversity (7) (Nei and 
Li 1979), net (Da) and absolute (Dxy) divergences (Nei 
1987), and to generate the haplotypes set including in- 
variant sites and gaps. POPART (Leigh and Bryant 2015) 
was used to construct a haplotype network using a medi- 
an-joining algorithm and default settings. 


Results 


Based on the comparisons carried out (detailed below), 
we confirmed that the examined specimens from the Pan- 
do stream (Figs 1, 2A, B), Santa Lucia system (Fig. 2C), 
and middle Rio Negro basin (Fig. 2E) correspond to new 
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records of D. pampeana, which extend its distribution to 
the southwest from the upper Negro basin (in straight-line 
distances: ~ 200 km to Middle Negro, ~ 280 km to Santa 
Lucia, and ~ 290 km to Pando) (Fig. 3). The list of exam- 
ined specimens of D. pampeana is presented in Table 1. 


Morphological comparisons 


The measurements of the examined specimens are sum- 
marized in Table 2. Comparing the morphometric data 
between the Pando group and the other groups, discrete 
differences between the ranges obtained were not detect- 
ed. Some tendencies based on the mean in some measure- 
ments, with partially overlapping ranges, were observed. 
The distance between the snout and dorsal-fin origin tend- 
ed to be slightly longer in the Pando, Santa Lucia and Mid- 
dle Negro groups than in the Upper Negro group (53.5- 
58.2% SL, mean = 55.7%+1.4 in Pando; 54.2—55.7% SL, 
mean 54.9%+1.1 in Santa Lucia; 54.2-57.6% SL, mean 
= 55.8%+1.0 in Middle Negro vs. 49.0-56.2% SL, mean 
= 52.7%+1.7 in Upper Negro). The Pando group present- 
ed a slightly smaller distance between the dorsal- and 
adipose-fin origins compared to the Upper Negro group 
(30.4-34.8% SL, mean = 33.0%+1.1 vs. 33.0-39.9% 
SL, mean = 37.1%+1.4), but almost similar to the other 
groups (33.6—-34.4 & SL, mean = 34.0+0.6 in Santa Lucia; 
32.0-36.1% SL, mean = 33.9%+1.0 in Middle Negro, re- 
spectively). In other measurements, as the caudal pedun- 
cle length and snout length, the Pando group tended to 
show greater mean values as follow (Pando, Santa Lucia, 
Middle Negro and Upper Negro, respectively): for cau- 
dal peduncle length 12.1-14.6% SL, mean = 13.6%+0.7; 
12.1-13.4% SL, mean = 12.8+0.9; 12.5—15.0% SL, mean 
= 13.4%+0.6; 8.9-13.4% SL, mean = 11.3%+1.0; and for 
snout length 19.3—22.3% HL, mean = 21.0%+0.8; 21.5— 
21.6% HL, mean = 21.5%+0.1; 20.1—22.1% HL, mean = 
20.8%+0.6; 16.3—22.1% HL, mean = 19.1%+1.6. 

Based on the consensus between the scree plot method 
and broken-stick model (Suppl. material 2), to ensure that 
did not discard biologically pertinent data, four eigenvec- 
tor elements were selected in the PCA on the size-correct- 
ed data, which accounted for 70.7% of the total variance 
(Suppl. material 3; Fig. 4A, B). Along the first axis 1n the 
PC1 vs. PC2 plot (Fig. 4A: explained 53.2% of the total 
variance), the Pando group was almost fully differentiat- 
ed from the Upper Negro group, but overlapped with the 
majority of the specimens of the Santa Lucia and Middle 
Negro groups. In the PC3 vs. PC4 plot (Fig. 4B explained 
17.6% of total variance), the groups appeared to overlap, 
and there was no clear distinction between them. PC1 
was most heavily loaded by the following measurements 
(Fig. 4A, Table 3): negatively by the snout to dorsal-fin 
origin (-0.5), snout to anal-fin origin (-0.3), caudal pe- 
duncle length (-0.3), snout to pelvic-fin origin (-0.2), and 
caudal peduncle depth (-0.2); and positively by the dor- 
sal- to adipose-fin origins (0.5) dorsal fin to caudal-fin 
base (0.2), dorsal-fin length (0.2), and anal-fin base length 
(0.2). PC2 was most influenced by positive variables such 
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Figure 1. Coloration in life of D. pampeana (A, B) from the Pando stream, Canelones Uruguay. Photo by J. Pfleiderer. 


Table 1. New records and material examined of D. pampeana. n= number of examined specimens. Group names corresponds with 
those described in the text. Accession numbers: OR533516* and OR533515**. 


Group 


Santa Lucia 
Santa Lucia 
Pando 
Pando 
Pando 
Pando 
Pando 
Middle Negro 
Upper Negro 


Upper Negro 
Upper Negro 
Upper Negro 
Upper Negro 
Upper Negro 


Upper Negro 


Upper Negro 


n 


10 


SL (mm) 


IB) 
23:9 
32.0=33.9 
30.6 
22:2=2019 
25.3-34.8 
252 
19.6-29.8 
29:9=33.6 


27.4-28.7 
24.3 
27.2-32.0 


27.3-29.1 
255 
25:9-33.6 


29.6 


Catalog 
number 


MHNM 1125 
MHNM 1189 
MHNM 812 
MLP 14443 
MLP 11444* 
MLP 11445 
UNMDP 5219** 
MHNM 4018 
UFRGS 8119 


UFRGS 8120 
UFRGS 8121 
UFRGS 8122 


UFRGS 8123 
UFRGS 8429 
UFRGS 8464 


UFRGS 28705 


Country 


Uruguay 
Uruguay 
Uruguay 
Uruguay 
Uruguay 
Uruguay 
Uruguay 
Uruguay 
Uruguay 


Uruguay 
Uruguay 
Uruguay 
Uruguay 
Brazil 


Brazil 


Brazil 


Locality 


Canelones, Rio Santa Lucia basin, Canelon Grande stream 
Canelones, Rio Santa Lucia basin, Canelon Grande stream 
Canelones, Canada de Ramos, Pando, Pando stream 
Canelones, Canada de Ramos, Pando, Pando stream 
Canelones, Canada de Ramos, Pando, Pando stream 
Canelones, Canada de Ramos, Pando, Pando stream 
Canelones, Canada de Ramos, Pando, Pando stream 
Durazno, marginal lagoon to Rio Yi, Estancias del Lago 
Cerro Largo, small stream at Route 26, ca. 59 km from 
Melo, between Sauce creek and Fraile Muerto creek 
Tacuarembé, Rio Tacuaremb6, at Route 26, Villa Ansina 
Rivera, Mazangano Bridge at Route 44 
Rivera, lateral puddles and Corrales creek, affluent of Rio 
Tacuaremb6, Route 27 
Tacuarembé, Caraguata creek, tributary to Rio 
Tacuarembé, Route 26, Las Toscas 
Rio Grande do Sul, Bagé, road between Acegua and 
Bagé, Rio Negro 
Rio Grande do Sul, Bagé, road between Acegua and 
Bagé, BR-153, Cinco Saltos creek, affluent of Rio Negro 
Rio Grande do Sul, Bagé, road between Acegua and 
Bagé, BR-153, Cinco Saltos creek, affluent of Rio Negro 


Latitude/Longitude 


34°29'14.00'S, 56°20'33.61"W 
34°29'14.70'S, 56°20'34.54"W 
34°43'39.01"S, 55°56'39"W 
34°44'19.2'S, 55°56'27"W 
34°42'12'S, 55°56'42.6'W 
34°44'19.2'S, 55°56'27"W 
34°42'12'S, 55°56'42.6'W 
33°21'47.16'S, 56°35'23.43'W 
32°17'39'S, 54°44'59"W 


JIE DO.33,5,.00- Lolo w, 
32°06'33'S, 54°40'08.6"W 
o1723. 20-9, 00 VOL4a W 
D2n09 29-3; 00-0 27a, 
31°28'37'S, 54°08'20"W 
31°36'53'S, 54°08'42"W 


31°36'53'S, 54°08'42"W 


Remarks 


1 c&s: 28.7 mm SL 


15 fully measured 


3 c&s: 30.4— 
31.5 mm SL 


holotype 
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Figure 2. Extern morphology of studied specimens of Diapoma pampeana. A. MLP 11443, male, 30.6 mm SL, Uruguay, Pando 
Stream; B. MLP 11445, female, 35.1 mm SL, Uruguay, Pando Stream; C. MHNM 1125, female, 31.3 mm SL, Uruguay, Canelon 
Grande Stream; D. MHNM 4018, male, 29.8 mm SL, Uruguay, marginal lagoon to Rio Yi; E. MHNM 4018, female, 29.3 mm SL, Uru- 
guay, marginal lagoon to Rio Yi. Photographs of the specimens from the Upper Negro are available in Ito et a/. (2022). Scale bar: 1 mm. 
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Figure 3. Geographic distribution of Diapoma pampeana. A. View within South America; B. View within Brazil and Uruguay. All 
studied specimens from the Pando (Circle), Santa Lucia (triangle), Yi (diamond), and Upper Negro (star; holotype represented by 
not-filled pattern) are depicted. Other records presented by Ito et al. (2022) are also included. 


Table 2. Comparative morphometric data obtained in the specimens studied of D. pampeana. SD = standard deviation; n = number 
of examined specimens; CI = confidence interval. Group names corresponds with those described in the text. 


Pando (n = 17) Santa Lucia (n = 2) Middle Negro (n = 15) Upper Negro (n = 35) 

Range MeanszSD CI95% Range MeansSD CI95% Range Mean+SD CI95% Range MeansSD CI95% 
Standard Length (mm) 22.2-35.1 28.9+3.8 27.1; 30.6 28.9-31.3 30.1+1.7 28.9; 31.3 23.3-29.8 26.641.7 25.8; 27.4 25.1-33.6 29.44+2.2 28.7; 30.1 
Percents of SL (%) 
Depth at dorsalfin 29.5-34.7 31.741.4 31.1; 32.4 30.1-32.2 31.141.5 30.1; 32.2 29.6-33.4 31.5+1.3 30.9; 32.1 27.7-33.7 30.941.2 30.5; 31.3 
origin 
Snout to dorsal-fin 53.5-58.2 55.7+1.4 55.1;56.4 54.2-55.7 54.941.1 54.2;55.7 54.2-57.6 55.8+1.0 55.4;56.3 49.0-56.2 52.741.7 52.3; 53.3 
origin 
Snout to pectoraHfin  24.7-27.3  26.2+0.9 25.8; 26.6 25.2-26.2 25.7+0.7 25.2; 26.2 25-28.3 26.4+0.8 26.0; 26.8 23.2-27.5 25.1+0.9 24.8; 25.4 
origin 
Snout to pelvic-fin 44.1-47.8 46.5+1.0 46.1; 47.0 44.7-45.0 44.9+0.2 44.7; 45.0 44.3-48.2 45.8+1.0 45.3:46.3 41.9-47 44.9+1.2 44.5; 45.3 
origin 
Snout to analfin origin 56.6-62.1 59.9+1.5 59.3; 60.7 57.9-58.5 58.2+0.4 57.9;58.5 56.6-60.6 58.7+1.2 58.2;59.3 53.8-60.6 58.3+1.5 57.8; 58.8 
Distance between 30.4-34.8 33.0+1.1 32.5; 33.5 33.6-34.4 34.0+0.6 33.6; 34.4 32.0-36.1 33.9+1.0 33.4; 34.4 33.0-39.9 37.1+1.4 36.6; 37.5 
dorsal- and adipose-fin 
origins 
Dorsal fin to caudal-fin 46.4-49.0 47.7+0.9 47.3;48.1 45.5-48.6 47.14+2.2 45.5;48.6 46.1-49.7 48.2+1.3 47.6; 48.9 45.9-54.8 49.04+1.9 48.2; 49.3 
base 
Dorsalfin length 20.7-24.8 23.0+0.9 22.6; 23.5 23.2-24.0 23.6+40.6 23.2; 24.0 22.7-25.1 23.940.7 23.6; 24.3 22.7-27.0 24.641.2 24.2; 25.0 
Dorsalfin base length 9.7-12.8 11.1+0.7 10.7;11.4 10.4-12.2 11.34+1.3 10.4;12.2 10.5-11.9 11.0+0.4 10.8;11.3 10.1-14.0 11.7+0.9 11.4; 12.0 
Pectoral-fin length 21.6-25.0 23.141.0 22.6; 23.6 22.2-23.2 22.740.7 22.2; 23.2 20.9-24.4 22.6+0.9 22.2; 23.1 21.9-25.8 23.840.9 23.5; 24.1 
Pelvic-fin length 11.1-14.7. 1340.9 12.5;13.4 12.2-12.3 12.2+0.1 12.2;12.3 11.2-13.7 12.640.8 12.2;13.0 12.0-15.7 13.9+0.9 13.5; 14.2 
Anal-fin base length 30.2-37.0 33.5+1.8 32.7; 34.3 31.4-35.1 33.3+2.6 31.4; 35.1 31.7-34.3 33.0+0.9 32.5; 33.4 32.1-36.6 34.5+1.3 34.0; 34.9 
Caudal peduncle depth 8.6-11.6 10.1+0.7 9.8;10.4 10.5-10.7 10.6+0.2 10.5;10.7 9.6-11.0 10.4+0.4 10.2;10.7 7.8-10.0 9.0+0.5 8.8;9.2 


Caudal peduncle 12.1-14.6 13.640.7 13.3;13.9 12.1-13.4 12.8+0.9 12.1;13.4 12.5-15.0 13.4+0.6 13.1;13.7 8.9-13.4 11.3+1.0 11.0; 11.6 
length 

Head length 21.9-25.4 23.7+40.9 23.3; 24.1 21.7-22.9 22.34+0.9 21.7; 22.9 22.5-25.3 23.54+0.8 23.1; 23.9 21.3-25.2 23.140.9 22.8; 23.3 
Percents of HL (%) 

Snout length 19.3-22.3 21.0+0.8 20.7; 21.4 21.5-21.6 21.5+0.1 21.5; 21.6 20.1-22.1 20.8+0.6 20.5; 21.1 16.3-22.1 19.1+1.6 18.6; 19.7 
Horizontal eye length 38.8-45.0 42.1+1.6 41.4; 42.9 41.8-42.8 42.3+0.7 41.8;42.8 40.2-44.5 42.3+1.3 41.7;42.9 39.5-46.8 43.5+1.6 43.0; 44.1 
Postorbital head 36.7-42.8 39.2+2.1 38.2; 40.1 38.3-40.2 39.3+1.3 38.3; 40.2 35.4-40.0 38.0+1.3 37.3; 38.6 35.6-43.1 39.24+1.9 38.5; 39.8 
length 

Least interorbital width 30.1-36.0 33.341.8 32.5; 34.1 34.2-34.9 34.5+0.5 34.2; 34.9 28.8-35.9 33.44+1.8 32.6; 34.3 28.6-38.7 32.24+1.9 31.6; 32.8 
Upper jaw length 32.8-44.5 36.54+2.9 35.1; 37.8 39.1-39.9 39.5+0.6 39.1; 39.9 32.7-37.1 35.341.4 34.7; 36.1 32.4-39.1 35.7+1.8 35.1; 36.3 
Dentary length 37.1-46.4 40.3+2.4 39.1; 41.3 39.1-39.8 39.4+0.5 39.1; 39.8 37.7-42.8 40.3+1.6 39.5; 41.0 36.5-42.2 38.94+1.6 38.3; 39.4 
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Figure 4. Most discriminant axes obtained from the PCA analyses performed using morphometric and meristic data of studied spec- 
imens of Diapoma pampeana (in each plot, the loadings are scaled to 90% of the PC scores). Size-corrected measurements: A. PC1 
vs. PC2 plot; B. PC3 vs. PC4 plot. Meristic data; C. PC1 vs. PC2 plot; D. PC3 vs. PC4 plot. Only these variables that most loaded the 
components are indicated as follows: E- depth at dorsal-fin origin; F- snout to dorsal-fin origin; G- snout to pelvic-fin origin; H- Snout 
to anal-fin origin; I- distance between dorsal- and adipose-fin origins; J- dorsal fin to caudal-fin base; K- anal-fin base length; L- 
caudal peduncle length; M- longitudinal scales; N- lateral-line scales; P- scales between lateral line-dorsal origin; Q- scales between 
lateral line-pelvic origin; R- circumpeduncular scales; S- predorsal scales; T- number of branched anal-fin rays; U- gill rakers on up- 
per limb of branchial arch; V- gill rakers on lower limb of branchial arch; W- number of maxillary teeth; X- number of dentary teeth. 


as the snout to anal-fin origin (0.3), depth at dorsal-fin 
origin (0.2), snout to pelvic-fin origin (0.2), and distance 
between the dorsal- and adipose-fin origins (0.2). PC3 
was most affected by the dorsal fin to caudal-fin base 
(-0.3), whereas PC4 was most strongly loaded by the 
anal-fin base (0.3). 

The cluster analysis showed that the four groups ana- 
lyzed were distributed into two large clusters (most boot- 
strap values were below 50). All the specimens of the Up- 
per Negro group (except three) were almost completely 
separated from the Pando, Santa Lucia, and Middle Ne- 
gro groups. In contrast, the specimens of the Pando group 
were not clustered separately, but were instead mixed 
mainly with the specimens of the Santa Lucia and Middle 
Negro groups (Suppl. material 4). 

The comparative results obtained in the meristic data 
are presented in Table 4. The first four components (ex- 
plained 64.4% of the total variance, Suppl. material 3) 
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were chosen as significant to analyze the variation in the 
meristic data, following the same criteria used for the 
morphometric data (Fig. 4C, D; Suppl. material 2). The 
Pando group was slightly differentiated from the Upper 
Negro group along the horizontal axis in the PC1 vs. PC2 
plot (accounted for 41.3% of the total variance, Fig. 4C), 
but overlapped almost completely with the other groups. 
In the PC3 vs. PC4 plot (explained 23.0% of the total 
variance, Fig. 4D), the Pando group was not separately 
distributed from the other groups along the axes. PC1 was 
most strongly influenced by the number of longitudinal 
scales (0.7) and predorsal scales (0.7) (Fig. 4C, Table 3). 
PC2 was strongly loaded by the number of dentary (0.7) 
and maxillary (0.5) teeth, and number of gill rakers on 
the upper limb of the first branchial arch (0.5) (Table 3). 
PC3 was mainly influenced by the number of pored lat- 
eral-line scales (-0.7), whereas PC4 was greatly affected 
by the number of branched anal-fin rays (-0.7). Tukey 
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Table 3. Loadings obtained from the PCA analyses using mor- 
phometric and meristic data. Percentages of variance are reported. 


Variables Components 
1 2 3 4 
Morphometric data: 41.5% 11.6% 9.2% 8.4% 
Depth at dorsal-fin origin -0.1 0.2 0.2 0.2 
Snout to dorsalfin origin 0.5 0.1 0.1 0.0 
Snout to pectoral-fin origin -0.1 0.1 0.0 0.0 
Snout to pelvic-fin origin 0.2 0.2 0.0 0.0 
Snout to analfin origin -0.3 0.3 0.1 0.1 
Distance between dorsal- and adipose-fin origins 0.5 0.2 0.1 0.2 
Dorsal fin to caudal-fin base 0.2 0.1 0.3 0.2 
Dorsalfin length 0.2 0.0 0.1 0.1 
Dorsal-fin base length 0.1 0.0 0.1 0.0 
Pectoral-fin length 0.1 0.1 0.0 0.0 
Pelvic-fin length 0.1 0.1 0.1 0.1 
Analfin base length 0.2 0.1 0.0 0.3 
Caudal peduncle depth 0.2 0.0 -0.1 0.0 
Caudal peduncle length 0.3 0.0 0.2 0.0 
Head length -~0.1 0.1 0.0 0.0 
Snout length -~0.1 0.0 0.0 0.0 
Horizontal eye length 0.0 0.0 0.0 0.0 
Postorbital head length 0.0 0.1 0.0 0.0 
Least interorbital width 0.1 0.0 0.0 0.0 
Upper jaw length 0.0 0.1 0.0 0.0 
Lower jaw length 0.1 0.1 0.0 0.0 
Meristic data: 27.0% 14.4% 13.2% 9.8% 
Longitudinal scales 0.7 -0.4 0.1 0.1 
Lateralline scales 0.3 0.0 0.7 0.0 
Scales between lateral line-dorsal origin -0.6 -0.4 0.2 -0.4 
Scales between lateral line-pelvic origin -0.6 0.2 0.2 0.5 
Circumpeduncular scales 0.5 0.3 0.6 0.4 
Predorsal scales 0.7 0.3 0.1 0.0 
Branched anal-fin rays 0.2 0.2 -0.4 0.7 
Gill rakers upper limb of branchial arch 0.5 0.5 -0.3 0.0 
Gill rakers lower limb of branchial arch 0.6 0.1 0.1 -0.2 
Maxillary teeth -0.4 0.5 -~0.1 -0.1 
Dentary teeth 0.0 0.7 0.5 0.1 


box plots of counts that most affected PCA and were most 
distinctive for the Pando group are presented in Suppl. 
material 5. The multivariate analyses performed on the 
morphometric and meristic data converged in coincident 
results that, despite having overlap between some indi- 
viduals, showed the population from the Pando stream to 
be somewhat distinctive morphologically from the speci- 
mens from the upper Rio Negro basin. The number of ver- 
tebrae of the Pando group was observed within the range 
of variation of the Upper Negro group (34 vs. 34-35). 


7/ 


The pigmentation pattern observed in the Pando group 
(Figs 1, 2A, B) was similar to that found in the Upper Ne- 
gro groups (Ito et al. 2022: figs 1-3), mainly character- 
ized by the vertically enlarged humeral spot, the narrow 
and conspicuous black line along the horizontal septum 
of body (in some specimens of the Pando group, it was 
observed to be somewhat silvery), the longitudinal black 
stripe extending posteriorly on the middle caudal-fin rays, 
and the presence of a small black blotch, restricted to the 
base of the middle caudal-fin rays. The specimens of the 
Middle Negro group (Fig. 2D, E) were observed to be sim- 
ilarly pigmented as the other groups, except for the cau- 
dal-fin blotch, which was not completely extended along 
the middle rays in some specimens. The examined speci- 
mens of the Santa Lucia group were found slightly faded 
(Fig. 2C), but some pigmentation characters of D. pam- 
peana as those aforementioned were found to be present. 


Molecular comparisons 


The genetic variation seen, based on the Tajima-Nei dis- 
tance between the Pando specimens and the specimens 
from the upper Negro basin of D. pampeana were found 
to be very low or nearly zero (< 0.002), even when boot- 
strapped (10000). In the NJ topology (Fig. 5), the two 
specimens of the Pando stream analyzed were placed to- 
gether with the remaining specimens of D. pampeana from 
the upper Negro basin and, particularly, as closely related 
to one specimen of that basin than between each other 
(Fig. 5). The Bayesian topology showed that D. pampeana 
was more related to D. guarani and D. obi within a com- 
mon clade with D. potamohadros and D. tipiaia (Fig. 6). 
In both the Bayesian and NJ trees, the analyzed specimens 
of Diapoma from the Pando stream were strongly placed 
together in the same clade with the specimens of D. pam- 
peana. The general pattern of interrelationships among the 
Diapoma species was found to be almost similar between 
both methods. The uncorrected pairwise mean distance 
obtained for D. pampeana ranged from 3.5 to 7.6% (Sup- 
pl. material 6: the lowest value with D. obi and the high- 
est value with D. itaimbe). The intraspecific variation of 
D. pampeana was observed to be varying from 0 to 0.2%. 
Only one specimen of the Upper Negro group (UFRGS 


Table 4. Comparative meristic data obtained for the studied specimens of D. pampeana. SD = standard deviation; n = number of 
examined specimens. Mean and mode values are reported. Group names corresponds with those described in the text. 


Pando 

Range Mean/ n Range 

Mo Mode+SD 
Longitudinal scales 35-38 =36.9/3740.9 17 35-37 
LateraHine scales 7-9 7.8/7+0.8 17 5-8 
Scales between lateral line-dorsal origin 5-5 5.0/5+0.0 17 5-5 
Scales between lateral line-pelvic origin 4-5 4.1/4+0.2 17 5-5 
Circumpeduncular scales 14-15 14.1/14+0.3 15 15 
Predorsal scales 12-15 13.6/1340.8 17 14-15 
Branched anal-fin rays 21-26 8=23.4/2341.4 17 22-26 
Gill rakers upper limb of branchial arch 7-10 8.3/8+0.8 14 7 
Gill rakers lower limb of branchial arch 14-18 14.9/1441.2 14 15 
Maxillary teeth 1-3 1.6/2+0.6 15 2 
Dentary teeth 7-11 8.8/9+1.1 15 12 


Santa Lucia Middle Negro Upper Negro 
Mean/ n Range Mean/ n Range Mean/ n 
Mode+SD Mode+SD Mode+SD 
36.0/N/A#1.4 2 35-39 37.0/3641.1 15 32-37 35.1/3641.2 35 
6.5/N/A#2.1 2 5-8 6.7/7+1.0 15 5-9 7.3/8+1.0 35 
5.0/N/At0.0 2 5-5 5.0/5+0.0 15 56 5.6/6+0.5 35 
5.0/N/A40.0 2 4-5 4.1/4+0.4 15 45 4.5/4+0.5 34 
N/A 1 14-15 14.3/1440.5 15 11-15 13.0/134+0.9 35 
14.5/N/A+0.7. 2 11-13 12.3/12+0.6 15 10-14 12.3/12+0.8 35 
24.0/N/A#2.8 2 21-25 22.8/2341.1 15 21-25 9 22.8/22+1.1 35 
N/A 1 7-8 7.5/7+0.5 15 6-9 7.1/7+0.9 34 
N/A 1 14-15 14.3/1440.5 15 13-15 13.9/14+0.8 35 
N/A 1 1-2 1.3/1+0.5 15 1-4 2.1/3+0.9 34 
N/A 1 7-11 9.0/9+1.0 15 6-14 9.0/9+1.3 35 
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Figure 5. Neighbor-Joining topology of analyzed Diapoma specimens based on Tamura-Nei model and COI sequence data. Boot- 
strap values (10000 replicates) are shown below the branches. SBL = 0.783. 
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Diapoma itaimbe UFRGS12717_TEC1453 
0.8834 Diapoma itaimbe UFRGS12084_TEC1231h 
Diapoma itaimbe UFRGS12634_TEC1369e 
Diapoma itaimbe isolate UFRGS12534_TEC1236e 
4 Diapoma itainmbe UFRGS12585_TEC1261h 
Diapoma itaimbe UFRGS12600_TEC1300e 
Diapoma itaimbe UFRGS12620_TEC1355g 
0.9996 Diapoma itaimbe UFRGS12651_TEC1383 
Diapoma itaimbe UFRGS16517_TEC2845c 
Diapoma itaimbe UFRGS12625_TEC1360h 
09282'__ Diapoma itaimbe UFRGS12627_TEC1362h 
Diapoma itaimbe UFRGS16567_TEC2899b 
Diapoma dicropotamicum UFRGS12727_TEC1465a 
0.999 Diapoma dicropotamicum UFRGS12727_TEC1465c 
Diapoma dicropotamicum UFRGS12727_TEC1465d 


0.9597 
1 
0.9483 


0.9626] Diapoma dicropotamicum UFRGS12727_TEC1465e 
Diapoma alburnum MCP21291 
areas Diapoma alburnum UFRGS12471_TEC1279 


7 Diapoma alburnum UFRGS10753_TEC1390 
07718 Diapoma alburnum UFRGS12593_TEC1269 
Diapoma alburnum UFRGS15029_TEC1810c 


Figure 6. Bayesian topology of phylogenetic relationships among the analyzed Diapoma species (comparing specimens of 
D. pampeana from the Pando stream and Upper Negro basin) based on COI sequence data. Numbers at nodes correspond to poste- 


rior probabilities. 


12642: TEC1377e) showed the greater p-distance (0.2%) 
in each comparison with the other specimens analyzed of 
D. pampeana (Suppl. material 6). 

In the polymorphism analysis comparing the Pando and 
Upper Negro specimens, 728 sites were analyzed (381: 
invariable; 347: with gaps or missing data), resulting in 
one polymorphic site (singleton) and two haplotypes (Hd 
= 0.286; variance = 0.034; standard deviation = 0.196). 
In general, the nucleotide diversity was extremely low (1 
= 0.00095; 0 = 0.00107; k = 0.286) for the samples ana- 
lyzed. No variation was found within the Pando samples 
(x = 0.00000; k = 0.000). For the Upper Negro samples, 
only one polymorphic site (monomorphic in the Pando 
samples) was found and, thus, their diversity was slightly 
greater (1 = 0.00105; k = 0.400). There were no observed 
shared mutations between these two populations. Regard- 
ing the divergence between the populations compared, the 
values obtained were low (Dxy = 0.00052; Da = 0.00000). 
The haplotype network showed a simple structure of two 
groups without well-defined geographic structure and in 
which one of them was mixed (Suppl. material 7). 


Comparative examined material 


Diapomaalburnum:. UFRGS 13309, 11, 33.4-56.0 mm SL. 
Diapoma guarani: MHNG 2366.99, holotype, 31.7 mm 
SL. Diapoma lepiclastum: MACN-ict 9682, 47, 29.3— 
42.0 mm SL. Diapoma obi: MLP 11312, 3, 29.5—35.6 mm 
SL. MACN-ict 9560, holotype, 52.6 mm SL. Diapoma 
uruguayense, MACN-ict 9681, 7, 31.6—-34.6 mm SL. 


Discussion 


The stevardiine species D. pampeana was recently de- 
scribed from several localities along the Rio Negro basin 
in Brazil and Uruguay (Ito et al. 2022). The freshwater 
fish fauna from the Pando stream is mainly known from 
studies focused on estuarine-influenced coastal waters 
(Plavan et al. 2010; Gurdek and Acufia-Plavan 2017). Al- 
though the body coloration was similar between the spec- 
imens of the Pando and Upper Negro groups, the former 
showed a striking morphometric and meristic divergence 
from the latter. Additionally, the intraspecific variation 
between the Pando and Upper Negro populations be- 
came much more subtle when compared with specimens 
from geographically intermediate areas such as the San- 
ta Lucia and Middle Negro basins (i.e. the specimens of 
these groups were slightly more similar to each other in 
the morphometric and meristic data than to the specimens 
of the Upper Negro). Frequently, in morphological com- 
parisons using multivariate methods, the populations ex- 
hibiting major differences in body shape correspond with 
those morphotypes that are located at the farthest distance 
geographically from the others (Lazzarotto et al. 2017; 
Vanegas-Rios et al. 2019; Rodrigues-Oliveira et al. 2023). 
In consequence, the analyzed specimens of D. pampea- 
na may be responding to a gradual pattern of divergence 
associated with spatial segregation, often observed in 
widespread species (Lazzarotto et al. 2017; Arroyave et 
al. 2019; Vanegas-Rios et al. 2019; Rodrigues-Oliveira et 
al. 2023). 
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In the morphological comparisons performed herein, 
the Pando group was greatly differentiated in the size-cor- 
rected PCA from the Upper Negro group along PC1, 
which was mainly influenced by the following distances: 
snout to dorsal-fin origin, snout to anal-fin origin, dor- 
sal- and adipose-fin origins, and caudal peduncle length 
(Fig. 4A). Some measurements such as the snout to dor- 
sal-fin origin, snout to anal-fin origin, snout to pelvic-fin 
origin, distance between the dorsal- and pectoral-fin or- 
igins, head length, snout length, and eye diameter have 
been found to be taxonomically informative to discrimi- 
nate among Diapoma species (Malabarba and Weitzman 
2003; Vanegas-Rios et al. 2018; Ito et al. 2022). Based 
on the meristic data, these groups were also found to be 
slightly divergent from each other along PC1. For this 
component, some main counts defining the variability 
were the number of longitudinal scales and predorsal 
scales (Table 3, Fig. 4C). Usually, the range is used as 
the main indicator to define the limits of the morpholog- 
ical variation among species in measurements or counts 
(Garavello et al. 1992; Aguirre et al. 2016; Lazzarotto et 
al. 2017; Arroyave et al. 2019; Vanegas-Rios et al. 2019; 
Malabarba et al. 2021). The main dilemma for defining 
these limits appears when data have varied degrees of 
overlapping between populations of study. 

Further statistical procedures are used as a complement 
to test if diverging tendencies in morphometric data are, 
or are not, significant (Lazzarotto et al. 2017; Arroyave et 
al. 2019; Vanegas-Rios et al. 2019; Rodrigues-Oliveira et 
al. 2023). For instance, it is frequently found that subtle or 
moderate intraspecific differences in body shape result in 
being statistically significant, but it may be also associat- 
ed with the degree of sensitivity involved in pairwise tests 
(e.g. means). Recently, the statistical potential behind 
morphometrically divergent patterns has been used to 
propose subspecies within gymnotids (Craig et al. 2017). 
Such criterion is not commonly used in modern ichthy- 
ology, so that the erection of infraspecific categories can 
be considered unnecessary to understand the intraspecific 
variation (e.g. clines) (Kottelat 1998; Kullander 1999). 
In the case of the Pando group, it seems unjustified to 
propose a new infraspecific category under the light of 
the current evidence described from morphological data, 
even more when it is geographically isolated and its phe- 
notypic variation is described and contextualized in the 
present contribution. In a similar case, geographically iso- 
lated coastal river populations of D. itaimbe that showed 
a Statistically significant difference in overlapping ranges 
of anal-fin ray counts were treated as structured isolated 
populations instead of separate species (Malabarba and 
Weitzman 2003; Hirschmann et al. 2015). 

When the discriminative tendencies are striking, as oc- 
curred here between the Pando and Upper Negro groups, 
including additional independent evidence, such as DNA 
data, can help to support the conclusion. The COI marker 
has played an important role in resolving taxonomic ques- 
tions in freshwater fishes (Pereira et al. 2013). In general, 
comparative studies using COI and morphology are deal- 
ing with cryptic species, species complexes or populations 
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with slight morphological variations (under a context of 
geographic isolation) (Serrano et al. 2019; Garavello et al. 
2021; Guimaraes et al. 2021; Malabarba et al. 2021; Aguil- 
era et al. 2022). The phylogenetic signal of the COI marker 
within Diapoma has been studied and used to propose in- 
traspecific and interspecific limits, as well as new species 
(Casciotta et al. 2012; Hirschmann et al. 2015; Ito et al. 
2022). However, the use of this marker by itself for recog- 
nizing species (e.g. as single locus without morphological 
or cytogenetic support) has its own methodological lim- 
its, and it is especially important to have this into account 
when dealing with complex groups of species or popu- 
lations (Castro Paz et al. 2014; Garcia-Melo et al. 2019; 
Klimov et al. 2019; Silva-Santos et al. 2023). This may 
imply that in some cases there is no guarantee of complete 
interspecific delimitation (e.g. Castro Paz et al. 2014). 

The phylogenetic comparison performed using the COI 
marker of all known species of Diapoma (except D. nandi 
from the Parana basin) recovered the Pando group as part 
of D. pampeana. It also demonstrates that the recognition 
of the Pando group as separate would make D. pampeana 
paraphyletic (Fig. 5: NJ topology). The variability of the 
p-distances calculated for the Pando group were observed 
to be lower than the average congeneric values (0O-0.2% 
vs. 1.3-8.0%), as it has been often reported in other 
characids (Pereira et al. 2011b; Garcia-Melo et al. 2019; 
Silva-Santos et al. 2023). Furthermore, the mean genetic 
distances were 3.5% and 4.3% between D. pampeana and 
its closest related species (D. obi and D. guarani, respec- 
tively). The intraspecific values obtained for D. pampea- 
na are within the conspecific variation reported in spe- 
cies of other characid groups such as Astvanax (from the 
Rio Paraguac¢u basin, mean = 0-1.7%); Hyphessobrycon 
Durbin, 1908 (from the Amazon basin, mean = 0-8.9%), 
and some stevardiines (Bryconamericus, Eretmobrycon 
Fink, 1976, Hemibrycon Gunther, 1864, Knodus Eigen- 
mann, 1911, and Piabina: 0-1.9%) (Pereira et al. 2011b; 
Castro Paz et al. 2014; Garcia-Melo et al. 2019; Sil- 
va-Santos et al. 2023). In barcoding studies, 2% threshold 
limit (at least 10 times the average conspecific values) has 
been used as the cutoff divergence value for delimiting 
species or molecular operational taxonomic units (Hebert 
et al. 2004; Ward 2009). However, an alternative thresh- 
old value of 1% has been considered for studying species 
complexes (Hubert et al. 2008; Pereira et al. 2011a; Sil- 
va-Santos et al. 2023). Therefore, these threshold values 
should be cautiously evaluated for each case in particular. 

The results obtained from the exploratory analysis using 
haplotype, allowed us to detect the potential presence of 
two haplotypes that were separated by a single mutational 
change. Additionally, the Pando group presented the same 
haplotype as the Upper Negro group, which again reinforc- 
es the great resemblance between both groups. Additional- 
ly, no well-defined lineages were detected in the molecular 
comparisons. However, this needs to be further investigat- 
ed so that within Diapoma, for instance, it has been found 
that D. itaimbe forms populations with well-defined struc- 
tural lineages associated with a coastal biographic pattern 
(Hirschmann et al. 2015). The recognition of species as 
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separately evolving metapopulation lineages is a unifying 
concept in defining species (De Queiroz 2007), and so far, 
we have no support from molecular data to separate the 
Pando group from D. pampeana. 

Freshwater fishes have limited their ability to disperse 
across brackish, marine or terrestrial barriers, being bio- 
logically restricted to water bodies after their formation, 
and thus, the disjunct geographic range associated with 
Species and populations across several basins might be ex- 
plained by river captures or dispersal favored by temporary 
connections (Albert and Reis 2011; Thomaz et al. 2017; 
Camelier et al. 2018; Cassemiro et al. 2023). These poten- 
tial explanations would be plausible for testing in D. pam- 
peana if its distribution was really widespread along most 
coastal drainages in Uruguay, as our findings suggest. 


Conclusion 


We concluded that the specimens from the Pando stream, 
despite the morphological divergence observed, can be 
classified as D. pampeana. We supported our decision 
based on the following arguments: 1) the deviations found 
on the morphometric and meristic data (e.g. PCA) are not 
enough to erect a new species and, as consequence, the 
intraspecific variability is increased; 2) the specimens of 
the Pando group were similarly pigmented as the speci- 
mens of the Upper Negro group (sharing the same diag- 
nostic pattern on the humeral mark, midlateral stripe, and 
caudal-fin pigmentation); 3) the COI-based phylogenetic 
procedures supported the placement of the Pando group 
within the genetic variation of the Upper Negro group 
(typical distribution of D. pampeana), and 4) based on 
the genetic distances, the Pando group was found to be 
genetically similar to the Upper Negro group, with p-dis- 
tances (O—0.2%) being lower than the mean distances ob- 
tained between each congener (1.3—8.0%). Additionally, 
the present work also confirmed the presence of D. pam- 
peana in the Yi (Middle Negro basin) and Santa Lucia 
river basins, based on morphological evidence. Although 
it was not possible to separate species, our results provide 
new information that can be further appreciated. For in- 
stance, it has been proposed that diverging populations 
can represent separate evolutionarily significant units, 
which should be conserved (Moritz 1994; de J. May-Itza 
et al. 2012; Stockwell et al. 2013; Berger et al. 2018). The 
geographic range of D. pampeana seems to be incom- 
pletely understood and might be more widely represented 
along the Rio Negro basin and other coastal drainages 
in Uruguay (Fig. 3, Suppl. material 8). Future studies 
may bring new insights into the population variation of 
the species and other phylogeographic patterns if more 
specimens and new localities are analyzed. The Pando 
population of D. pampeana constitutes a remarkable ex- 
ample of an isolated population that is morphologically 
divergent (in the morphometric and meristic data) from 
the geographically most distant conspecific population 
(upper Rio Negro basin), but that shares a high degree of 
genetic resemblance with it. 
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